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Abstract 

It is shown that for a given bipartite density matrix and by choosing a suitable separable set 
(instead of product set) on the separable-entangled boundary, optimal Lewenstein-Sanpera (L-S) 
decomposition can be obtained via optimization for a generic entangled density matrix. Based on 
this, We obtain optimal L-S decomposition for some bipartite systems such as 2(8)2 and 2(8)3 Bell 
decomposable states, generic two qubit state in Wootters basis, iso-concurrence decomposable 
states, states obtained from BD states via one parameter and three parameters local operations 
and classical communications (LOCC), Werner and isotropic states, and a one parameter 

3 <8 3 state. We also obtain the optimal decomposition for multi partite isotropic state. It is 
shown that in all 2 (8 2 systems considered here the average concurrence of the decomposition 
is equal to the concurrence. We also show that for some 2 (8 3 Bell decomposable states the 
average concurrence of the decomposition is equal to the lower bound of the concurrence of 
state presented recently in [Buchleitner et al, quant-ph/0302144], so an exact expression for 
concurrence of these states is obtained. It is also shown that for d <8 d isotropic state where 
decomposition leads to a separable and an entangled pure state, the average I-concurrence of 
the decomposition is equal to the I-concurrence of the state. 

Keywords: Quantum entanglement, Optimal Lewenstein-Sanpera decomposi- 
tion, Concurrence, Bell decomposable states, LOCC 

PACS Index: 03.65.Ud 
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1 Introduction 

In the past decade quantum entanglement has been attracted much attention in connection with 
theory of quantum information and computation. This is because of potential resource that entan- 
glement provides for quantum communication and information processing [1, 2, 3]. By definition, 
a bipartite mixed state p is said to be entangled if it can not be expressed as 

P = S Wi ^ ^ 2) ' Wi - °' J2 Wi = 1 ' 

i i 

where pp and pf^ denote density matrices of subsystems 1 and 2, respectively. Otherwise the 
state is separable. 

The central tasks of quantum information theory is to characterize and quantify entangled states. 
A first attempt in characterization of entangled states has been made by Peres and Horodecki et 
al. [4, 5]. Peres showed that a necessary condition for separability of a bipartite system is that 
its partial transpose be positive. Horodecki et al. have shown that this condition is sufficient for 
separability of composite systems only for dimensions 2 2 and 2 <g) 3. 

Having a well justified measure to quantify entanglement, particularly for mixed states of a 
bipartite system, is indeed worth, and a number of measures have been proposed [3, 6, 7, 8]. 
Among them the entanglement of formation has more importance, since it intends to quantify the 
resources needed to create a given entangled state. 

Another interesting description of entanglement is Lewenstein-Sanpera (L-S) decomposition 
[9, 10]. Lewenstein and Sanpera have shown that any bipartite density matrix can be represented 
optimally as a sum of a separable state and an entangled state. They have also shown that for two 
qubit systems the decomposition reduces to a mixture of a mixed separable state and an entangled 
pure state, thus all entanglement content of the state is concentrated in the pure entangled state. 
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This leads to an unambiguous measure of entanglement for any two qubit state as entanglement 
of pure state multiplied by the weight of pure part in the decomposition. The strategy of Refs. 
[9, 10] is based on the fact that for a given set V = {|e a , f a )} of product states belonging to the 
range of density matrix p, one can subtract separable density matrix p* = J2 a A a -P Q (not necessary 
normalized) with A a > such that bp = p — p* > 0. 

In Ref. [9], the best separable approximation (BSA) has been obtained numerically in case of 
two qubit Werner state by choosing a set of several hundred P a -projectors. Some analytical results 
is also obtained for special states of two qubit states [11]. Further, in [12] BSA of a two qubit 
state has been obtained algebraically. They have also shown that in some cases the weight of the 
entangled part in the decomposition is equal to the concurrence of the state. In Ref. [13] we have 
obtained optimal L-S decomposition for a generic two qubit density matrix by using Wootters basis. 
It is shown that the average concurrence of the decomposition is equal to the concurrence of the 
state. 

In this paper we obtain optimal L-S decomposition for some bipartite systems. Here we obtain 
optimal decomposition for a given density matrix p by choosing suitable separable set S in which 
p s G S. This approach is different from the others in the sense that optimal decomposition is 
obtained for a given separable set S instead of product set V. Also this approach is geometrically 
intuitive as it will be explained in section 4 by providing a bunch of interesting bipartite systems 
such as, 2 (g) 2 and 2 ® 3 Bell decomposable states, a generic two qubit state in Wootters basis, 
iso-concurrence decomposable states, states differing from BD states via one parameter and three 
parameters local operations and classical communications (LOCC), d (g> d Werner and isotropic 
states, and a one parameter 3 3 state. We also provide the optimal decomposition for multi 
partite isotropic system. As a byproduct we show that in all 2 2 systems considered here the 
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average concurrence of the decomposition is equal to the concurrence. We also show that for some 
2 ® 3 Bell decomposable states for which entangled part of the decomposition is only a pure state 
the average concurrence of the decomposition is equal to the lower bound of the concurrence of 
state presented recently in Ref. [14], consequently an exact expression for concurrence of these 
states is given. In the case of d ® d isotropic state we show that the average I-concurrence of the 
decomposition is equal to the I-concurrence of the state. 

The paper is organized as follows. In section 2 we, briefly, review concurrence as presented in 
[8]. In section 3 we first review Lewenstein-Sanpera decomposition for bipartite density matrix, 
then a new prescription for finding optimal decomposition is presented. Some important bipartite 
examples is considered in section 4. The paper is ended with a brief conclusion in section 5. 

2 Concurrence 

In this section we review concurrence of two qubit mixed states as introduced in [8] . The generalized 
concurrence is also reviewed, briefly. 

2.1 Wootters's Concurrence 

From the various measures proposed to quantify entanglement, the entanglement of formation has 
a special position which in fact intends to quantify the resources needed to create a given entangled 
state [3]. Wootters in [8] has shown that for a two qubit system entanglement of formation of a 
mixed state p can be defined as 




(2-1) 
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where H(x) = — xlnx — (1 — x) In (1 — x) is binary entropy and concurrence C{p) is denned by 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }, (2-2) 

where the A« are the non-negative eigenvalues, in decreasing order, of the Hermitian matrix R = 
\J\fPP\fP where the spin-flipped state p is defined by 

P=(<J y <g> (Ty)p*((7y <g) (Ty), (2-3) 

where p* is the complex conjugate of p in a standard basis such as {|00) , 1 01) , 1 10) , |11)} and cr y 
represent Pauli matrix in local basis {|0) , |1)}. 

Consider a generic two qubit density matrix p with its subnormalized orthogonal eigenvectors 
\vi), i.e. p = J2i \ v i) ( v i\- There always exist a decomposition [8] 

p = J2\xi)( Xi \ (2-4) 

% 

where Wootters's basis \xj) are defined by 

4 

l^> = E^>*>' for * = 1,2,3,4, (2-5) 
j 

such that 

( Xi | xj) = {UTU T )ij = XiSij, (2-6) 

where r^- = (vi \ Vj) is a symmetric but not necessarily Hermitian matrix. The states |xj), which 
are going to be used in our notation, is defined as 

\x' i ) = ^, for i = 1,2,3,4. (2-7) 



2.2 I-concurrence 

Several attempts to generalize the notion of concurrence for arbitrary bipartite quantum system 
have been made already [15, 16, 17]. Among them the so-called I-concurrence [17] is defined in 
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terms of universal-inverter superoperator which is a natural generalization to higher dimensions of 
two qubit spin flip. I-concurrence of a joint pure state of a d,A <8> ds system is defined by Rungta 
et al. [17] 

C(m = y/2(l - tr(pl)) = ^2(l-tr(p|)), (2-8) 
where pa = tred^*) (V*!) an d Pb is defined similarly 

3 Lewenstein-Sanpera decomposition 

According to Lewenstein-Sanpera decomposition [9], any bipartite density matrix p can be written 

as 

p = Xp s + (1 - X)p e , A €[0,1], (3-9) 

where p s is a separable density matrix and p e is an entangled state. The Lewenstein-Sanpera (L-S) 
decomposition of a given density matrix p is not unique and, in general, there is a continuum set 
of L-S decomposition to choose from. However, Lewenstein and Sanpera in [9, 10] have shown that 
the optimal decomposition is unique for which A is maximal. Furthermore they have demonstrated 
that in the case of two qubit systems p e reduces to a single pure state. 

The idea of Refs. [9, 10] is based on the method of subtracting projections on product vectors 
from a given state, that is, for a given density matrix p and any set V = {\e a , f a )} of product 
states belonging to the range of p, one can subtract separable density matrix p* = J2 a h-aPa (not 
necessary normalized) with all A a > such that dp = p — p* > 0. Separable state p* provides 
the optimal separable approximation (OSA) in the sense that trace Tr(p* s ) < 1 is maximal and 
entangled part p e is called edge state, a state with no product vectors in the range [18]. Lewenstein 
and Sanpera provide the conditions that trace Tr(p*) is maximal. 
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In this paper we will deal with L-S decomposition from different point of view. Our approach 
is based on the fact that the set of separable density matrices is convex and compact [19, 20]. This 
follows from the fact that any separable density matrix p s G S can be written as a finite convex 
combination of pure product states. The set of all Hermitian operators acting on the Hilbert 
space Hi <8> Ti.2 constitute a Hilbert space (called Hilbert-Shcmidt space) with a real inner product 
< A, B >= Tr(A^B). The set of density matrices p are defined as Hermitian positive semi-definite 
and the trace one matrices form subset V of H-S space which is compact and convex [20]. Let V 
denotes the set of all pure product states. V is tensor product of two spheres which are compact in 
the finite dimensional case. So V is also compact [19]. The set of all finite convex combinations of 
product states V is defined as the convex hull of V, i.e. S = coirvV, and convex hull of a compact 
set V is also compact, so the set of separable density matrices is compact [19]. 

Based on the above fact we obtain optimal L-S decomposition for some bipartite systems. For 
a given density matrix p, we choose a suitable separable set S C S on the separable-entangled 
boundary, and express p as a convex combination of separable state p s G S and an arbitrary 
entangled state p e , i.e. p = Xp s + (1 — X)p e . Then we evaluate A and provide the conditions that A 
is maximal under the restrictions that p s is in the separable set S and maintaining the positivity of 
the difference p — Xp s , i.e. p e remains nonnegative. To this aim we allow p s to move on the surface 
defined by S, and simultaneously search for the p e with corresponding maximal A. This restricts 
p e to some entangled states and gives p s as a function of p and restricted p e . The only matter that 
should be noticed in choosing the set S for which p s G S, is that all states on the line segment 
connecting p s and p, i.e. p e = ep s + (1 — e)p for < e < 1, must be entangled. This guarantees 
that thus obtained decomposition is indeed maximal. In all examples considered in this paper we 
will see that the rank of p e is less than rank of p. This means that p e is an edge state with no 
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product vectors in its range as pointed out in Ref. [18]. Moreover in the case of two qubit system 
it is shown that p e reduces to pure entangled state as we expect from the results of Refs. [9, 10]. 
For these systems p s is defined as a function of p and concurrence of entangled pure state. To make 
our consideration more clear, we provide some examples in the next section. 

4 Some important examples 

In this section we obtain optimal decomposition for some categories of states, namely, 2® 2 Bell de- 
composable (BD) states, a generic two qubit state in Wootters basis, iso-concurrence decomposable 
states, some 2 ® 2 states obtaining from BD states via one parameter and three parameters LOCC 
operations, 2 03 Bell decomposable states, d®d Werner and isotropic states, a one parameter 3(g> 3 
state and finally multi partite isotropic state. 

4.1 2 ® 2 Bell decomposable states 

We begin by considering the 2 ® 2 Bell decomposable (BD) states. A BD state acting on H 4 = 
H 2 ® H 2 Hilbert space is defined by 



4 



4 




< pi < 1, 




(4-10) 



i=i 



i=i 



where \tpi) are Bell states given by 



IV'i) = k + ) = ^ ! (|oo) + |ii)), 
IV>3> = |^ + ) = -Ul 01 > + l 10 >)> 



iv> 2 > = \r) = ^(ioo) 

IV>4> = |YO = -Ul 01 > 



in)) 



|io». 



(4-11) 



(4-12) 



A BD state is separable iff pi < ^ for all i 



= 1, 2, 3, 4 [21]. In the following we consider the case that p 



is entangled for which p\> \. To obtain optimal L-S decomposition we choose p s = J2i=i Pi IV'i) (V^ 
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with p\ = i as boundary separable state and p e = J2i=iPi Itfi) (V'il- Inserting these equations into 
the decomposition given in Eq. (3-9) we get 

Pi = ApJ + (1 - A)j# for i = 1,2,3,4. (4-13) 

From Eq. (4-13) we get A = and aTF 7 = -§n > where C = 2 Pl - 1 and C" = 2p'{ - 1 are 

concurrence of p and p e , respectively. This means that in order to obtain optimal decomposition, 
i.e. having maximal A, we should require that C" takes its maximal value, where this happens as 
long as p'2 = p'i = p'l = 0, i.e. p e is pure entangled state. Considering the above arguments we get 
for A, p s and p e the following results 

A = l-C, Pe = |Y>i)(Vi|, 

(4-14) 

Pi = h Pj = T fori = 2, 3, 4. 
Equation (4-14) simply shows that average concurrence of the decomposition is equal to the con- 
currence of state, i.e. (1 — A)C(|^)) = C. 

4.2 A generic two qubit state in Wootters's basis 

In this subsection we obtain optimal L-S decomposition for a generic two qubit density matrix 
by using Wootters basis. In Ref. [13] we have shown that a generic two qubit density matrix 
p = J2i ^i\ x i) ( x i\ w ith corresponding set of positive numbers Aj and Wootters's basis \x'j) can 
be obtained from a Bell decomposable state with the same set of positive numbers Aj but with 
different Wootters's basis via SO(4, c) transformation. It is also shown that local unitary trans- 
formations correspond to 50(4, r) transformations, hence, p can be represented as coset space 
SO(4, c)/SO(4, r) together with positive numbers \. 

Now in order to obtain optimal L-S decomposition we choose p s = J2i K \ x 'i) ( x 'i\ with A'i — A' 2 — 
A3 — A4 = as boundary separable state and p e = J2i A'/ \ x 'i) ( x i\- Inserting these equations into 
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the decomposition given in Eq. (3-9) we get 

Ai = AA- + (1 - A)A" for i = 1,2,3,4. (4-15) 

From Eq. (4-15) we get A = C 'q7, C and = > where C = X\ - A 2 - A 3 - A 4 and 
C" = X" — X f 2 — A3 — A4 are concurrence of p and p e , respectively. This means that in order to 
obtain optimal decomposition, i.e. having maximal A, we should require that C" takes its maximal 
value, which happens as long as X% = A3 = A 4 ' = 0, i.e. p e is pure entangled state with concurrence 
A". Considering the above arguments we get for A, p s and p e the following results 

A = i-§> p e = X'i\x' 1 )(x' 1 \, 



(4-16) 

_ _ c>.- •) •) 1 

A ' j A 

Equation (4-16) simply shows that average concurrence of the decomposition is equal to the con- 



A / = a 2+ A3+a 4 ; wAi for j = 2,3,4. 



currence of state, i.e. (1 — X)C(\tp)) = C. The decomposition (4-16) is in agrement with results 
obtained in Ref. [13] 

4.3 Iso-concurrence decomposable states 

In this section we define iso-concurrence decomposable (ICD) states, then we give their separability 
condition and evaluate optimal decompsition. The iso-concurrence states are defined by 

= cos0|OO> +sin0|ll)), |0 2 > = sin0|OO) -cos0|ll», (4-17) 
\(j> 3 } = cos 0|O1) +sin0|lO)), |0 4 > =sin0|Ol) -cos0|lO)). (4-18) 

It is quite easy to see that the above states are orthogonal thus span the Hilbert space of 2 ® 2 
systems. Also by choosing = f the above states reduce to Bell states. Now we can define ICD 
states as 

4 4 
P = I>I&M&I' 0<pi<l, = (4-19) 



i=l i=l 
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These states form a four simplex (tetrahedral) with its vertices defined by p± = 1, P2 = 1, P3 = 1 
and p4 = 1, respectively. 

Peres-Horodeckis criterion [4, 5] for separability implies that the state given in Eq. (4-19) is 
separable if and only if the following inequalities are satisfied 



(Pi - V2) < \J 4p 3 P4 + (P3 ~ Pa) 2 sin 2 29, 



(P2 - Pi) < \J 4p 3 P4 + (P3 ~ Pa) 2 sin 2 29, 



(P3 ~Pa) < y 4pip 2 + (pi - P2) 2 sin 2 20, 
(P4 - P3) < y/ 4pip 2 + (Pi - P2) 2 sin 2 26. 



(4-20) 
(4-21) 
(4-22) 
(4-23) 



Inequalities (4-20) to (4-23) divide tetrahedral of density matrices to five regions. Central regions, 
defined by the above inequalities, form a deformed octahedral and are separable states. In four 
other regions one of the above inequality will not hold, therefor they represent entangled states. 
Bellow we consider entangled states corresponding to violation of inequality (4-20) i.e. the states 
which satisfy the following inequality 



(Pi - P2) > \J 4p3P4 + (P3 - P4) 2 sin 2 29. 



(4-24) 



All other ICD states can be obtain via local unitary transformations. Now we will obtain concur- 
rence of ICD states. Following the Wootters protocol given in subsection 2.1 we get for the state p 
given in Eq. (4-19) 

/ \ 

— pism29 y/PiP2 cos 29 

v / PiP2 cos 29 P2sin26> 

P3sin2# — y/p3P4 cos 29 

— y/P3P4 cos 29 — j»4sin20 



T = 



(4-25) 
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Now it is easy to evaluate Aj which yields 



Ai, 2 = \ (±(Pi ~ P2) sm 26 + yj A Pl p 2 + (pi - p 2 ) 2 sin 2 26^ , 
A 3 ,4 = 3 (±(pa - Pi) sin 26* + ^ 4p 3 p 4 + (P3 - P4) 2 sin 2 2dj . 



(4-26) 



Thus one can evaluate the concurrence of ICD states as 

C=(pi-p 2 ) sin 26 - \/ 4p 3 P4 + (P3~ Pa) 2 sm 2 26. (4-27) 

It is worth to note that thus obtained concurrence is equal to the amount of violation of inequality 
(4-24). Note that the concurrence of an ICD state can be written as 

A u - A 22 - ^J(A 33 + A44) 2 - 4^ 2 4 , (4-28) 

where Aij denote matrix representation of ICD state in Bell basis, that is 

A u = -fa +p 2 + (pi -p 2 )sin26»), A 22 = -( Pl + p 2 - ( Pl - p 2 ) sin 26), (4-29) 

-433 = 2 (P3 + Pi + (P3 - P4) sin 2(9) , ^44 = ^ (P3 + P4 - (P3 - P4) sin 20) , (4-30) 

Au = -(pi -p 2 )cos2#, ^34 = -(pa -p 4 )cos2fl. (4-31) 

Now in order to obtain optimal L-S decomposition we parameterize p s like ICD state with matrix 
elements A'- (in Bell basis) which are defined like A^ except for pi and 6 which are replaced with 
p\ and 6', respectively. We also choose p e similar to p with matrix elements A"j parameterized with 
p'( and 6". For simplicity the rank of p e is considered to be two, namely p'l = p'l = 0. Using these 
consideration together with Eq. (3-9) we get 

A ij = \A! ij + (l-\)A!l j , (4-32) 

Taking into account the fact that p s is boundary separable state with zero concurrence and using 
Eq. (4-28), we get A = C 'qv, C and ^4-, where C and C" are concurrence of p and p e , respectively. 
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Obviously we observe that A becomes maximal when p e is a pure entangled state. Considering this 
fact and setting p 2 ' = we arrive at 

Pi + P2 + (pi ~ P2) sin 29 = A (pi + p' 2 + (pi - p' 2 ) sin 29') + (1 - A)(l + sin 29"), (4-33) 

P1+P2- (Pi ~ P2) sin 29 = A (pi + j/ 2 - (pi - p' 2 ) sin 29') + (1 - A)(l - sin 29"), (4-34) 

P3+P4 + (Pa - P4) sin 26* = A (p' 3 + p 4 + (p' 3 - p 4 ) sin 20') , (4-35) 

P3 + P4 ~ (P3 - P4) sin 26* = A (p 3 + j/ 4 - (p'3 - p' 4 ) sin 29') , (4-36) 

(pi - P2) cos 2# = A(pi - p' 2 ) cos 20' + (1 - A) cos 29", (4-37) 

(P3 - Pa) cos 26> = A (p'3 - p' 4 ) cos 20'. (4-38) 

In order to solve above equations we consider two cases separability, 
i) case 1: 



First let us consider the case that 9,9' ^ j. In this case Eqs. (4-33) to (4-38) yield to 

9 = 9' = 9", 

A = 1 - (pi - p 2 ) sin 29 + ^/(p 3 +Pa)(P5+Pg), ( 4 ~ 39 ) 
pl = Elz § z Al j pj. = ^ for j = 2,3,4. 

This case corresponds to results of Ref . [22] . 
ii) case 2: 

Now let us consider the case that 9 = j, i.e. p is Bell decomposable state. The only nontrivial 
solution of Eq. (4-38) is p' 3 = p' A . Equations (4-35) and (4-36) show that this restricts the density 
matrix to p 3 = p±. Combining all, we arrive at the following p s for decomposition 

tan 29' = ?L±^zl tan 29", A = , (4-40) 
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1 / j_ C 1- Pl - P2 



/ / P3 

Ps=P4 = y- 



(4-41) 
(4-42) 



where C = 2p\ — 1 is concurrence of p. The separability of density matrix p s implies that p\ > 
for all i (recall that the separability condition has been already imposed over p s by putting its 
concurrence equal to zero) . So p\ should satisfy the following condition 



sin 20" > 



{pi+P2)C 
PiC + p 2 



(4-43) 



This condition also guarantees positivity of A. It is worth to emphasis that this case involves the 
result of Ref. [23] as a special case. There authors have obtained the optimal decomposition for a 
special kind of BD states, namely a specific Werner state with p\ = | (of course in their treatment 
they take singlet state ^4) as dominant pure state in Werner state, i.e. P 4 = |). 



4.4 One parameter LOCC operations 

A generic two qubit density matrix p can be represented in Bell basis as p = YAY^ where Y € 
SO(4, c)/SO(4, r) and A = diag(Ai, A2, A3, A4) [13]. Here we consider the case that Y is a one 
parameter matrix as 

/ \ 

cosh# isinh^ 



Y = 



— isinh# cosh# 



1 



1 



(4-44) 
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thus 



= YAyt = 



V 



(4-45) 



/ 



Ai cosh 2 9 + A 2 sinh 2 9 i(X 1 + A 2 ) sinh 9 cosh 9 
-«(Ai + A 2 ) sinh cosh Ai sinh 2 9 + A 2 cosh 2 9 
A 3 

A 4 

Obviously, normalization condition leads to (Ai + A 2 ) cosh 29 + A3 + A4 = 1. We choose p s in the 
same form as p, i.e. p s = Y'A'Y' j where A' = diag(A' 1 , A' 2 , A 3 , A' 4 ) and Y' is defined as Y but here 
is replaced with 9'. Now in order to obtain optimal L-S decomposition we have to get a generic 
density matrix for p e . After doing so, it can be easily seen that Eq. (3-9) requires that p e has also 
the same form as p and p s , i.e. p e = Y"A"Y' /J< where A" = diag(A' 1 / , A 2 , A3, A4) and Y" is defined 
as Y but with 9" instead of 9. Inserting the above equations in Eq. (3-9) we get 



yay^ = A(y / A / y ,t ) + (1 - A)(y"A"y ,/t ). 



(4-46) 



Now multiplying Eq. (4-46) by Y" T and Y"* , respectively from left and right and using the 
orthogonality of Y" we get 



(Y" T Y)A(Y^Y"*) = \{Y" T Y')A\Y' ] Y"*) + (1 - A) A", 



(4-47) 



where it can be written as 

(Ai cosh 2 {9 - 9") + A 2 sinh 2 (9 - 0")) = A (a' x cosh 2 (9' - 9") + A' 2 sinh 2 (9' - 6>")) + (1 - A)A' 1 / , 

(4-48) 

(Ai sinh 2 {9 - 9") + A 2 cosh 2 {9 - 0")) = A (a' x sinh 2 (9' - 9") + A' 2 cosh 2 (9' - 0")) + (1 - A)A 2 ', 

(4-49) 

A 3 = AA^ + (1 - A) A3, (4-50) 
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A 4 = XX'i + (1 - X)X'i, (4-51) 

(Ai + A 2 ) sinh 2(9 - 9") + X(X[ + A' 2 ) sinh 2(9' - 9") = (4-52) 

Subtracting Eqs. (4-49), (4-50) and (4-51) from Eq. (4-48) and using the fact that p s is boundary 
separable state, hence having zero concurrence, i.e. X'± — X' 2 — X' 3 — X' 4 = 0, we get A = C "qiF , -§^n = 
-^jjrz > where C = X± — A2 — A3 — A4 and C" = X'[ — X 2 — A3 — A4, are concurrence of p and 
p e , respectively. This shows that maximal A is achieved when X% = A3 = A4 = 0, i.e. p e is pure 
entangled state with concurrence X". Implying the above results we can solve Eqs. (4-48) to (4-52) 
for A and p s where we get 

A = l-Ccosh2fl", (4-53) 



tmh2( * 9 n, (A 1+ A 2 )sinh2(g-n 
tanh2(0 -9)- — - — _ — _ -, 



(4-54) 



V 1 / (A 1 + A 2 )cosh2(fl-Q-C \ 

Al ' 2 " 2A { cosh2(^-n ± (As + A4) ) ' (4 " 55) 

^ = y. for j = 3, 4, (4-56) 

where in Eq. (4-53) we have used X" = cos ^ 2 e" which follows from normalization condition of p e . 
Finally from the positivity conditions for A and A« we see that the following inequalities should hold 

c^<<<I, co sh2 (,-,"><^ + I ^. (4-57, 

Note that thus obtained decomposition is not a special case of the decomposition considered in 
subsection 4.2. There we considered the case that all p, p s and p e were expressed in same Wootters 
basis. Here their Wootters basis parameterized differently, namely 9, 9' and 9" respectively. The 
optimal decomposition given by Eqs. (4-53) to (4-56) involves some interesting special cases as 
follows 
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case i) 9 = 9': In this case from Eqs. (4-53) to (4-56) we get 9" = 6, which yields the results 
of subsection 4.2 for a one parameter Wootters basis. 

case ii) 9 = 0, 9" / 0: This case leads to optimal decomposition of a BD state in terms of the 
non maximal entangled pure state. This case also can be considered as generalization of the result 
ofRef. [23]. 

case iii) 9 / 0, 6" = : This case leads to the optimal decomposition of a one parameter 
LOCC transformed BD state in terms of maximal entangled pure state. 

4.5 Three parameters LOCC transformed BD states 

Now we consider the case that p can be obtained from BD states via three parameters LOCC 
transformation as p = YAY^ with [13] 

cosh 9 cosh £ cosh <\> + sinh 9 sinh i (cosh 9 cosh £ sinh <p + sinh 9 cosh (f>) i cosh 9 sinh £ ^ 



Y 



— i (sinh 9 cosh £ cosh <j> + cosh 9 sinh <f>) sinh 9 cosh £ sinh (f> + cosh 9 cosh <j) sinh 9 sinh £ 
— i sinh £ cosh sinh £ sinh <f> cosh£ 

1 

(4-58) 

where normalization condition leads to 



/ 



Tr(p) = ( (Ai cosh 2 4> + A 2 sinh 2 cosh 2 f + A 3 sinh 2 £ + (Ai sinh 2 4> + A 2 cosh 2 ) cosh 29 

+ ( Ai cosh 2 ^ + A 2 sinh sinh 2 £ + A 3 cosh 2 £ + ( Ai + A 2 ) cosh £ sinh 29 sinh 20 + A 4 = 1 . (4-59) 

We choose below p s in the same form as p, i.e. p s = Y'A'Y'* where A' = diag(A' 1 , A' 2 , A' 3 , A4) and 
Y 1 are defined as Y but here 6*, £ and are replaced with 9', £' and 0'. Now to obtain optimal 
L-S decomposition we should take a generic density matrix for p e . It can be easily seen that 
Eq. (3-9) requires that p e has also the same form as p and p s . So we get p e = Y"A"Y"* where 
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A" = diag(A' 1 / , A 2 , A3, A4) and Y" is denned as Y but here 9, £ and <f> are replaced with 9", £" and 
<f)" . By using the above considerations and Eq. (3-9) we get 

FAY 1 " = X(Y'A'Y' j ) + (1 - \){Y"M'Y" ] ). (4-60) 

Now multiplying Eq. (4-60) by Y" T and Y"*, respectively from left and right and using the 
orthogonality of Y" we get 

{Y" T Y)A(yW"*) = \(y" T Y')A'(Y*Y"*) + (1 - A) A", (4-61) 

Subtracting three last diagonal elements of matrix equation (4-61) from the first one and using the 
fact that p s has zero concurrence, i.e. A^ — A' 2 — A 3 — A 4 = 0, we get after some algebraic calculations 
A = C " C J, C and ^4- = > where C = X\ — A2— A3— A4 and C" = X'[— A 2 — A3— A4, are concurrence 
of p and p e , respectively. This shows that maximal A is achieved when A 2 = A 3 = A4 = 0, i.e. p e is 
pure entangled state with concurrence A'/. Considering the above results we can write Eq. (4-60) 
as 

p u = \ p ' n + (1 - A) X" (cosh 6» cosh £ cosh + sinh sinh 0) 2 , (4-62) 
P22 = Ap' 2 2 + (! - ^Mi (sinh 6> cosh £ cosh ^ + cosh 9 sinh <p) 2 , (4-63) 
P33 = Ap 33 + (1 - A)A? sinh 2 £ cosh 2 <f>, (4-64) 
pu = Ap 44 , (4-65) 
p 12 = Xp' 12 + (1 - A) A" ((cosh 2 £ cosh 2 <p + sinh 2 </») sinh 26 + cosh £ cosh 29 sinh 2</>) , (4-66) 
p 13 = Ap' 13 + (1 - A) A'/ (cosh 6> cosh 2 sinh 2£ + sinh sinh £ sinh 2^>) , (4-67) 
p 23 = Ap' 23 + ( 1 - A) A'/ (sinh 9 cosh 2 sinh 2£ + cosh 9 sinh £ sinh 2(f)) , (4-68) 
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where 

Pu = (Ai (cosh cosh £ cosh<^> + sinh sinh 0) 2 + A2 (cosh cosh £ sinh</> + sinh cosh 0) 2 

+A 3 (cosh sinh £) 2 ) , (4-69) 

P22 = (sinh 6* cosh £ cosh ^ + cosh sinh 0) 2 + A2 (sinh cosh £ sinh <fi + cosh cosh (f>) 2 

+A 3 (sinh sinh £) 2 ) , (4-70) 

p 33 = (Ai sinh 2 £ cosh 2 <p + A 2 sinh 2 £ sinh 2 + A 3 cosh 2 , (4-71) 

P44 = A 4 , (4-72) 

p 12 = ( ^Ai (cosh 2 £ cosh 2 ^ + sinh 2 <j^j + A 2 (cosh 2 £ sinh 2 + cosh 2 + A 3 sinh 2 sinh 20 

+ (Ai + A 2 ) cosh £ sinh 2(j) cosh 20) , (4-73) 

p 13 = cosh 2 ^ + A 2 sinh 2 </> + A 3 ) cosh 9 sinh 2£ + (Ai + A 2 ) sinh 9 sinh £ sinh 2^ , (4-74) 

p 23 = ((Ai cosh 2 + A 2 sinh 2 + A 3 ) sinh 9 sinh 2£ + (Ai + A 2 ) cosh 9 sinh £ sinh 20') , (4-75) 

and • are defined in same form as pij but here all parameters are expressed in terms of prime 
parameters. After tedious but straightforward calculations we arrive at the following results for p s 

tanh f - ~ F Sinh °' + G COSh °' 4 76) 

tan * ~ (p! +p2-A) sinh 20' + 1 cosh 20' ' ( " ' 

, -F cosh 0' + G sinh 0' 

tanh2£' = = 5 =- -, 4-77 

pi cosh 0' + p 2 sinh 9' + p% — \ {A cosh 20' — 2? sinh 29' + B + 2D) 

tanh2 ,/_ F sinh 9' - G cosh 0' 

tan ~ sinh £'(AA 3 + (pi + p 2 - A) cosh 20' - p 3 + E sinh 20' + D) ' ( " ' 

,/ 1 / —i 7 cosh 0' + G sinh 0' , 9 „. . , 9 „. „ 1.. . „„. ^ . , „ „^n\ 

X'= — [ Pl cosh 2 0' - p 2 sinh 2 0' + P 3 + - (A cosh 20' - E sinh 29' + B -2D)) 

2A V sinh2£' 2 / 

(4-79) 
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A; = ( ^2^ ( AA 3 + (Pi + ^2 - 4) cosh 20' -P 3 + E sinh 20' + £>) + AA' 3 + Pl - p 2 - p 3 - B + D 

(4-80) 

A' 2 = ( cos h 2( ^/ ( AA 3 + (Pi + P2 - ^4) cosh 29' -P 3 + E sinh 20' + D) - AX' 3 - Pl + p 2 + p 3 + B - D 

(4-81) 

^4 = y, (4-82) 

where 

A = (1 - A) A'/ ((cosh 2 £ cosh 2 </> + sinh 2 0) cosh 2(9 + cosh £ sinh 29 sinh 2^) , (4-83) 

B = (1 - A) A? (cosh 2 £ cosh 2 - sinh 2 </>) , (4-84) 

D = (1 - A) A'/ sinh 2 £ cosh 2 0, (4-85) 

£ = ( 1 - A) A" ( (cosh 2 £ cosh 2 + sinh 2 0) sinh 29 + cosh £ cosh 29 sinh 2^>) - p i2 , (4-86) 

F = (1 -A) A" (cosh 9 cosh 2 </> sinh 2£ + sinh 9 sinh f sinh 2</>) - pi 3 , (4-87) 

G = (1 - A) A'/ (sinh cosh 2 sinh 2£ + cosh 9 sinh f sinh 2</>) - p 23 . (4-88) 

The parameters 0' and £' are obtained by solving the Eqs. (4-76) and (4-77), using the remaining 
equations we can determine the parameters of p s in terms of parameters of p and p e . Note that 
the one parameter density matrix which was considered in previous subsection can be obtain from 
three parameters one by setting <p = <j>' = <p" = £ = £' = £" = 0. One can see that the equations 
in one parameter case is solvable and we can express the parameters of separable and entangled 
parts in L-S decomposition in terms of parameters of density matrix p which is the reason for its 
separated consideration in previous subsection. 
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4.6 2 <g) 3 Bell decomposable state 

In this subsection we obtain optimal L-S decomposition for Bell decomposable states of 2 (g) 3 
quantum systems. A Bell decomposable density matrix acting on 2 3 Hilbert space can be 
denned by 

6 6 

P = I>W*>M' 0< Pi <l, £> = 1, (4-89) 

i=i i=i 

where |Y>j) are Bell states in if 6 = ff 2 <g) ff 3 Hilbert space, defined by: 



hM = -±=(|ll> + |22», |^> = -^(|11>-|22» > 

|^ 3 ) = ^(|12) + |23», |V>4> = ^(|12) - |23», (4-90) 

|^5> = ^(|13> + |21» > |^ 6 ) = i=(|13)-|21)). 

It is quite easy to see that the above states are orthogonal and hence it can span the Hilbert space 
of 2 <g) 3 systems. From Peres-Horodeckis [4, 5] criterion for separability we deduce that the state 
given in Eq. (4-89) is separable if and only if the following inequalities are satisfied 

(pi - P2) 2 < (j>3 +Pi){Pb +P§), (4-91) 
(P3-P4) 2 < (P5+P6)(Pi +P2), (4-92) 

(P5 -P&? < (Pi +P2)(P3 +P4)- (4-93) 

In the sequel we always assume without loss of generality that p\ > P2, Pz > Pa and p^ > p§. 
Recently in Ref. [14] an analytical lower bound of concurrence of any 20 K mixed state is derived 
as 



C(p)> r£C*(pM), (4-94) 

V i>3 
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where p( l fi are unnormalized states restricted to 2 ® 2 subsystems under projection operators pto) 
as 

p (ij) = p(ij) pP {ij) i =I 2 ®(\i)(i\ + \j)(j\), (4-95) 

and C(pfe')) are Wootters concurrences of corresponding restricted 2 <g) 2 density matrices. For our 
2 Cg) 3 Bell decomposable state we get 

C(p {12) ) = max{0, pi - P2 - v / (P3 + P4)(P5 + P6)}, (4-96) 

C(/?( 23 )) = max{p 3 - p 4 - V / (Pi + P2)(P5 + P 6 )}, (4-97) 
C(/> (13) ) = max{ P5 - P6 - ^( Pl + p 2 )( P3 + P4 )}. (4-98) 

It is interesting to note that each Wootters concurrence given in Eqs. (4-96) to (4-98) corresponds 
to separability conditions given in Eqs. (4-91) to (4-93), respectively. Now in order to obtain 
optimal L-S decomposition for BD state given in Eq. (4-89) we choose p s = J^iPil^i) (V'il an d 
p e = J2iPi IV'i) (V'il- We also assume without loss of generality that p s lies on the separable- 
entangled boundary defined by (all other cases where p s lies on other surfaces can be treated 
similarly) 

P'l ~ p'i = ^(P's+P'M+Pe)- (4-99) 

Moreover p s must satisfies the other two separability conditions (4-92) and (4-93). This means that 
entangled state p violates separability condition (4-91), i.e. we have 

Pi > P2 + \J(P3 +Pi)(P5 + Pe). (4-100) 

However, two other inequalities (4-92) and (4-93) may be violated simultaneously. Taking into 
account the above considerations and Eq. (3-9) we get after some elementary calculations the 
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following equation 

(i - A) 2 ({pl-p'2) 2 - (p'I+p'IM+pD) 

-(1 - A) (2(pi - mWi-p'i) ~ (P3 +P4)(P5 +pl) - (P5 +Pe)(P3 +P4)) 

+ ((Pi -P2) 2 - (P3+P4)(P5+Pe)) =0. (4-101) 

Below in the rest of this subsection we will use Eq. (4-101) to calculate A for some possibile values 
of p'J, i = 1, 2, 6: as follows 
i) p" = 1: 

In this case Eq. (4-101) gives the following results 

A = l-pi-p2 + ^/(P3+ Pa)(P5 + Pe), Pe = |V»i) (V'll , (4-102) 



//, " P$ = ^ for j = 2, ■ • -,6. (1-L03) 



Furthermore p s must satisfies the separability conditions (4-92) and (4-93) which leads to the 
following restrictions for p 

{P3 - P4) 2 < (P5 + Pe) ( 2^2 + V(P3 +P4)(P5 + Pe)) 

v 7 (4-104) 

(P5 - Pe) 2 < (P3 +P4) (2p 2 + V(P3 +P4)(P5 + Pe)) 

By using Eq. (4-100) one can see that conditions (4-104) are stronger than separability conditions 
(4-92) and (4-93), that is in this case only separability condition (4-91) is violated by p. It is 
worth to mention that for these states we are enable to give exact expression for concurrence. 
As concurrence C(p) is defined as the infimum over all possible pure state decompositions, no 
decomposition can have average concurrence smaller than C{p). Since the decomposition given by 
Eqs. (4-102) and (4-103) constitute a maximal entangled pure state and a separable state p s , 
it follows that its average concurrence is equal to the weight of entangled part, namely (1 — A). On 
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the other hand for entangled states restricted by equations (4-100) and (4-104) we get C(p^) > 
and C(p^) = C(pW) = 0. This means that the lower bound is equal to (1 — A), i.e. 

C(p) = (1 - A) = pi - p 2 - ^J(P3+P4)(P5+P6)- (4-105) 

ii) P'l+P2 = l- 

In this case by performing optimization procedure Jp- = Jpr = in Eq. (4-101) (under 
constraint p'[ + p' 2 ' = 1), we can see that thus obtained equations from optimization procedure 
restrict the density matrix p to rank four one, namely p^ = p4 = or ps = p$ = 0. Under this 
circumstances we get A = C q,F and ^4z, where C and C" are concurrence of p and p e , respectively. 
This means that maximum A happens when p' 2 ' = which reduces to results of previous case. 

iii) p'[ + pg = 1 : 

After optimization procedure with the constraint p'[ + p 3 ' = 1 we get 

A = 1 - (pi - p 2 ) - (p3 + Pi) ~ |(P5 + Pe) 

(4-106) 

„/ _ 2p 2 -p 5 -p 6 / _ p 5 +P6-4p 4 / _ Pi f ■ _ 9 4 r fi 

where the following inequalities should be imposed in order p s to be separable state 

2 {pa ~ \{Pb + Pe)) < (j>5 + Pe) {p2 - \(P5 + Pe)) , 

2(P5 "Pe) 2 < (P5 +P6) (P2 " |(P5 +P6)) , ( 4 " 107 ) 
4^4 < P5 +P6 < 2p 2 - 

iv) pj+pg = 1: 

Analogue to the case p'[ + p 3 ' = 1 we get 

A = 1 - (pi -p 2 ) - j(p 3 +Pa) - (P5 + Pe), 

Pi — 2A ' ™3 — 4A ' Pj A ' IOT °* i t % °> 



(4-108) 
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whit restrictions 

2(P3 - Pi) 2 < (P3 + P4) (p2 - i(P3 + Pa)) , 

2 (pe - g(P3 + p 4 )) 2 < (P3 + P4) (j»2 - |(P3 + P4)) , ( 4 " 109 ) 
4j>6 < P3 +P4 < 2p 2 - 

v) pi'+pg + i^l: 

In this case it follows from optimization that rank p should be four, namely p^ = p% = 0. Under 
this conditions we get 

A = 2p 2 , 

(4-110) 

Pi = P2 = h P3 =P4= P5 = P6 = °- 

4.7 Werner states 

The Werner states are the only states that are invariant under U ®U operations. For d®d systems 
the Werner states are defined by [24] 

Pf = d^d ((d " f)I + (df ~ 1)F) ' _1 - f ~ (4 " ni) 

where / stands for identity operator and F = lu) 0'*|- It is shown that Werner state is 
separable iff < / < 1. Now to obtain L-S decomposition for Werner states we choose Pf=o as 
separable part and pf as entangled state , i.e. pf = Xpf=o + (1 — \)pf- Then from Eq. (3-9) we 
get A = ^-jr~ and ^ = -j^ < 0, that is A is maximum when /' = —1. Using the above results we 
get 

A = / + 1 ' P ' = W^T) {I - F) - ( " 12) 



Optimal L-S decomposition for some bipartite systems 27 

4.8 Isotropic states 

The isotropic states are the only ones that are invariant under U ®U* operations, where * denotes 
complex conjugation. The isotropic states of d <8> d systems are defined by [25] 

p F = ^I-(l-\^+)(^+\) + F\ij + ){^ + \, 0<F<1, (4-113) 

where = 4^ J2i I") is maximally entangled state. It is shown that isotropic state is separable 
when < F < ^ [25]. Now in order to obtain optimal L-S decomposition we choose boundary 
isotropic separable state with F = 1/d as separable part and ppr as entangled state where we get 

^ = ^dF^-i an< ^ W 7 = d (dF~-l$ — ^' * na * * s ' ^ ^ s max i mum when F' = 1. Using the above results 
we get 

A = ^TT^ p e =|^ + )(V + |. (4-114) 

It is interesting to stress that the average I-concurrence of the decomposition (4-114) is equal to 
the I-concurrence of the state obtained in Ref. [26]. By using Eq. (2-8) one can easily see that 
C(\ip + }) = \/2(l — 1/d), which can be used to evaluate average I-concurrence of the decomposition 

(l-A)C(|Y^)) = y^(V-i), for i<F<l, (4-115) 

which is equal to the I-concurrence of isotropic states which has been obtained in [26]. 

4.9 One parameter 3 ® 3 state 

Finally let us consider a one parameter state acting on H 9 = H 3 ® H 3 Hilbert space as [27] 

2 ol 5 — ol 

p* = ^ <V> + | + + — 2 < a < 5, (4-116) 
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where 

|^ + ) = ^ (|11) + |22) + |33)), 

a + = i(|12)(12||23)(23| + |31)(31|), (4-H7) 

<7_ = i(|21> (21||32) (32| + |13) (13|) . 
p a is separable iff 2 < a < 3, it is bound entangled iff 3 < a < 4 and it is distillable entangled state 

iff 4 < a < 5 [27]. To obtain L-S decomposition for p a we choose boundary separable state with 

a = 3 as p s and p e = p a /. After some calculations we get A = and ^ = ( 3 "~^ > 0. So the 

optimal L-S decomposition is achieved by choosing a' = 5 and we get 

A = ^, p e = 2 -\^)(^ + \ + \a + . (4-118) 
4.10 Multi partite isotropic states 

In this subsection we obtain optimal L-S decomposition for a n-partite d-levels system. Let us 
consider the following mixture of completely random state po = I/d n and maximally entangled 
state |V' + ) 

p( S ) = (l- S )-^ + s |V> + )<V + |, 0< S <1, (4-119) 

where / denotes identity operator in cP-dimensional Hilbert space and = Ylf=i !«*■■■*)■ 
The separability properties of the state (4-119) is considered in Ref. [28]. It is shown that the 
above state is separable iff s = sq = (1 + d n ~ l ) 1 . 

Now to obtain optimal L-S decomposition we choose p(so) & s separable part and p(s') as en- 

'— HA s— (l+d n_1 ) _1 

tangled part. By using Eq. (3-9) we get A = 8/ _ (1 ^. dw -i ) -i and ^ = ( a /_(i +d n-i ) -i) :j ' This means 
that the maximum A achieved when s' = 1, so we get 

A= (i-.Ki + < r-) | ftHW| . (4 . 120) 
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5 Conclusion 

We have shown that for a given bipartite density matrix and by choosing a suitable separable set on 
the separable-entangled boundary, optimal Lewenstein-Sanpera decomposition can be obtained via 
optimization over a generic entangled density matrix. Based on this , optimal L-S decomposition is 
obtained for some bipartite systems. We have obtained optimal decomposition for some bipartite 
states such as 2 2 and 2 ® 3 Bell decomposable states, generic two qubit state in Wootters basis, 
iso-concurrence decomposable states, states obtained from BD states via one parameter and three 
parameters LOCC operations, d® d Werner and isotropic states, a one parameter 3 ® 3 state and 
multi partite isotropic state. It is shown that in all 2 <g) 2 systems considered here the average 
concurrence of the decomposition is equal to the concurrence. We also obtain exact expression for 
concurrence of some 2(g) 3 BD states. In the case of d®d isotropic states it is shown that the average 
I-concurrence of the decomposition is equal to the I-concurrence of the states. We conjecture that 
for all optimal decomposition that entangled part is only a pure state, the average I-concurrence 
of the decomposition is equal to the I-concurrence of the state. 
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